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ANALISE MATEMÁTICA IV 

FICHA SUPLEMENTAR 4 

EQUAÇÕES DIFERENCIAIS LINEARES 



Formas canónicas de Jordan 



(1) 



Para cada uma das matrizes A seguintes, determine uma forma canónica de Jordan 


J, e uma matriz de mudança de base S tal que A = SJS' 1 . 


(a) (b) 


A = 


" -1 1 " 
-1 


A = 


' 1 -2 " 

2 1 




(c) (d) 


A = 


3 1 " 
-1 1 


A = 


' 1 0" 
1 1 
1 -1 2 





Resolução: 

(a) A matriz A é um bloco de Jordan, portanto J = A e 



S = 



1 
1 



(b) Os valores próprios da matriz são as soluções de 

det(A-XI) = ^^ (l-A) 2 + 4 = ^> A = 1 ± 2% . 
Conclui-se que uma forma canónica de Jordan de A é 



J = 



1 + 2» 

1 — 2i 



Os vectores próprios associados a 1 + 2i são os vectores 



que verificam 



-2t -2 

2 -2i 



-ia = b . 



Uma base do espaço próprio de 1 + 2i é constituída pelo vector 



í'i 



Os vectores próprios associados a 1 — 2% são os vectores 



que verificam 



2í -2 
2 2i 



ia = b . 



Uma base do espaço próprio de 1 — 2i é constituída pelo vector 



t'2 



Portanto uma matriz de mudança de base que põe A em forma canónica é 

1 1 



S = 



-i i 
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(c) Os valores próprios da matriz são as soluções da equação 



A 1 

-1 1-A 



= O 



X 2 - 4A + 4 = O 



(A - 2) 2 = O 



Logo A tem apenas um valor próprio, A = 2. Uma vez que a matriz não é igual a 21 
(onde I é a matriz identidade), o espaço próprio tem necessariamente dimensão 1 e 
portanto a forma canónica de Jordan de A é 



J = 



2 1 
O 2 



Os vectores próprios de A são os que verificam (A — 2I)v = O, ou seja 



Uma base dos vectores próprios é formada, por exemplo, pelo vector 

1 
-1 

Um vector próprio generalizado é um vector w que satisfaz (A — 2I)w = v. Podemos 
tomar, por exemplo, 

" 1 



te 



O 



Conclui-se que uma matriz de mudança de base que põe A em forma canónica é 

S = 



1 1 
-1 O 



(d) Os valores próprios da matriz são as soluções da equação 



1-A 
1 
1 



O 
-A 
-1 2 



-A(l-A)(2-A) + (1-A) = O ^^ 

(1-A)(-A(2-A) + 1) = O ^ 

(1_A)(1-A) 2 = O <=> 

(1-A) 3 = 0. 

Logo A tem apenas um valor próprio, A = 1. Os vectores próprios são os que verificam 



c = b — a 



Uma base do espaço próprio é constituída pelos vectores 









1 


-1 1 


1 


-1 1 





a 




' " 




b 


= 







c 








Vi = 



V2 = 



que se obtiveram fazendo a = l,b = 1 e a = —1, b = O, respectivamente. Conclui-se 
que uma forma canónica de Jordan de A é 



J = 



1 O O 
O 1 1 
O O 1 



AMIV - FICHA SUPLEMENTAR 4 3 

Falta achar um vector próprio generalizado correspondente a terceira coluna de J . 
Esse vector será uma solução de 

(A - I)w = v 

onde v é um vector próprio correspondente ao valor próprio 1 que gera o espaço das 
colunas da matriz (A — I). Uma vez que 

" O 

A-I= 1-11 

1 -1 1 

o espaço das colunas é gerado por v\ + v%- Assim, procura-se uma solução de 

(A — I)w = Vl + V2 ■ 

Toma-se, por exemplo, 



w = 



Portanto, uma matriz de mudança de base que põe A em forma canónica é 

"10 1 

S= 1 1 

1 o 

onde a primeira coluna é o vector próprio v\, a segunda coluna é o vector próprio 
v\ + V2 e a terceira coluna é o vector próprio generalizado w associado a v\ + v%- 

D 



(2) 



Para cada uma das matrizes A seguintes, determine uma forma canónica 
J, e uma matriz de mudança de base S tal que A = SJS' 1 . 


de Jordan 


(a) (b) 




A = 


2-22" 
1 
1 


A = 


2 
-1 
1 



-6 

4 




(c) (d) 




A = 


"20 1 " 
3-1 
1 1 


A = 


2 1 

2 

-1 


1 



3 




(e) (f) 




A = 


111" 

2 
-1 1 3 


A = 


" -2 4 
-2 4 

-4 4 


2 1 


4 
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Exponencial de Matrizes 



(3) 



Para cada uma das matrizes A seguintes, 


determine e At . 


(a) 


(b) 




A = 


' -TT " 

2tt 






A = 


"11" 
1 1 




(c) 


(d) 




A = 


" 1 2 ] 
1 






A = 


"4 5 " 
-2 -2 




(e) 








5 1 







A = 


5 


1 









5 













Resolução: 

(a) A exponencial de uma matriz diagonal, At, é a matriz diagonal cujas entradas são as 
usuais exponenciais escalares das entradas correspondentes em At. Deste modo, 

s -tí o 



At 







,2ttí 



(b) Os valores próprios de A são dados por 

(1-A) 2 -1 = ^^ A 2 -2A = 

^^ A = ou A = 2 . 

Os vectores próprios associados ao valor próprio são dados por 
Os vectores próprios associados ao valor próprio 2 são dados por 



11 
1 1 




a 
b 


= 


' " 




-1 1 
1 -1 



Portanto, 



,,At 



1 1 
-1 1 



e OÍ 

n „2t 



a = b . 



1 _i 
f I 

2 2 



S 



2 

-l + e 2i 



-1 + e 2 
lie 2 



(c) A matriz A tem apenas um valor próprio, A = 1, o qiva/ íem um espaço próprio de 
dimensão 1. Um vector próprio é uma solução de (A — I)v = 0, isto é, 

2 







v = 



Por exemplo pode-se tomar 



v = 



AMIV - FICHA SUPLEMENTAR 4 5 

Um vector próprio generalizado w, obtém-se resolvendo a equação (A — I)w = v, 
isto é, 

' O 2 



O O 



w = v 



Por exemplo, pode-se tomar 



w = 



Relativamente, a base (v,w), a transformação linear representada por A é dada por 
um bloco de Jordan, J, para o valor próprio 1, ou seja, 



A = 



"10" 




"li" 




10 


L \_ 




1 




2 



Logo, a exponencial é 



e At = S 



é té 
e* 



S~ x = 



s- 1 

ê 2te t 
e* 



Comentário: Notando que \A é um bloco de Jordan para o valor próprio \, podia- 



se, em alternativa, calcular 



Í4í 

e2 M 



e 2 te 2 



e^ 



e avaliar em 2í. 



<> 



(d) Os valores próprios de A são dados por 

(4-A)(-2-A) + 10 = <= 



A 2 -2A + 2 = 
A= l±i . 



Os vectores próprios (complexos) associados ao valor próprio A = 1 + i são dados por 



3-í 5 
-2 -3-i 



a 
ò 



o = — - — a . 
5 



Os vectores próprios para o valor próprio 1 — i são dados por 



3 + i 5 
-2 -3 + t 



Portanto, 



-,At 



1 1 

í— 3 — 3— i 



,(l+i)í Q 



o é 1 -^ 



h ~ 3 ~ i 

o = — - — a . 



1-3» 

2 


-5»" " 

2 


1+3» 
2 


5» 
2 - 



(e it +e- it )-3i(e it -e- it ) 
2 

i(e lt — e~ lt ) 



-5i(e H -e- u ) 
2 

(e it + e- it )+3i(e it -e- it ) 



= e 



cos t + 3 sin t 5 sin í 

—2 sin í cos í — 3 sin í 
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Comentário: Se X é valor próprio complexo da matriz real A com o vector próprio 
v, então o seu complexo conjugado também é valor próprio de A e o vector complexo 
conjugado de v é um vector próprio correspondente: Av = Xv => Av = Xv => 
Av = Xv. Esta observação permite simplificar cálculos. <0> 



(e) A matriz é um bloco de Jordan 3 por 3. 



e At = 



= 5í 



o 



te 



5í 
,,5í 



te 



D 5t 
òt 



.,5í 



D 



Sistemas de Equações Lineares 



(4) 



Considere a matriz 

A = 


"12" 
2 1 




(a) Quais são os valores próprios de Al 


(b) Quais são os vectores próprios de Al 


(c) Determine uma matriz de mudança de base, S, que diagonaliza A, e determine 


a sua inversa, S 1-1 . 


(d) Calcule e At . 


(e) Determine a solução do seguinte problema de valor inicial: 




m 
. 2/2 . 


= A 


S/i 

. 2/2 _ 


com 


" 3/i (0) n 

. y 2 (o) . 


= 


" 5 n 
7 




(f) Determine a solução geral da seguinte equação diferencial: 




S/i 

_ 2/2 . 


= A 


. 2/2 . 




(g) Escreva duas funções y : IR — > IR 2 que constituam uma base do espaço vectorial 


das soluções da equação da alínea anterior. 



Resolução: 

(a) Os valores próprios são os zeros do polinómio característico: 

p(X) = det(A - XI) = (1 - A) 2 - 4 . 

p(X) = ^> A 2 -2A-3 = ^> A=-louA = 3. 
Os valores próprios de A são -1 e 3. 

(b) Um vector v é vector próprio de A associado ao valor próprio X se e só se Av = Xv, 

ou seja, se e só se (A — XI)v = 0. Em componentes, escreve-se v = 
Equação para os vectores próprios associados ao valor próprio —1: 



-a = b . 



2 2" 
2 2 




a 
b 


= 


" " 




Os vectores próprios associados ao valor próprio A = — 1 são 



v = 



com 



a e 



AMIV - FICHA SUPLEMENTAR 4 
Equação para os vectores próprios associados ao valor próprio 3: 

^=> a = b . 



Os vectores próprios associados ao valor próprio A = 3 são 

aeR. 



v = 



a 
a 



com 



(c) Mudando para uma base de vectores próprios de A, a transformação linear fica dia- 
gonal. Tome-se, por exemplo, a matriz de mudança de base 



S = 



1 1 
-1 1 



cujas colunas são vectores próprios de A associados aos valores próprios -1 e 3. A 
mudança inversa é 



S~ l = 



í 

2 J 



Então tem-se 



A= S 
(d) De acordo com a alínea anterior, 



-1 O 
O 3 



.At 



= s 



O e 



5^ = 



s- 1 . 



2 

-t,„3t 



2 

-Í-LZ.3Í 



(e) A solução do problema de valor inicial dado é 



y(t) 



-.At 



2/1 (0) 

2/2(0) 



-*-Lz.3i 



-i_Lz,3i 



~ f + 6e 3í 

-t + g e 3í 



Vi G 



(f) /4 solução geral da equação diferencial dada é 

"* O 



y(í) = 5 



O e 



:.',/ 



Cl 
C2 



cie * + C2e 3í 
— cie - * + C2e 3í 



Vi g 



onde ci,C2 G R. 
(g) Compõe-se uma base para o espaço vectorial das soluções da equação da alínea 
anterior com as colunas da matriz 

-t 



S 



e - O 
O e 3í 



e"* e 3í 
-e"* e 3í 



ou se/a, coíT7 as funções x\ : 



Xl(t) = 



e x 2 



x 2 (t) = 



dadas por 

e 3t 

,;ít 
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De facto, x\(t) e X2(t) são funções linearmente independentes e qualquer solução 
y(t) da equação da alínea anterior é da forma y(t) = c\xi{t) + 02^2 (í) para algum 
cieKe algum c 2 G R. 

D 



(5) 



Considere a matriz 

A = 


"31" 

-1 1 






(a) Calcule e At . 




(b) Determine a solução do seguinte problema de valor inicial 






S/i 


= A 


S/i 

. 2/2 . 


com 


" 2/i(0) " 
„ 3/ 2 (0) 


= 


1 " 
-1 




(c) Determine a solução geral da seguinte equação diferencial 




í Vi = 3yi +y 2 +e 2t 




\ V2 = -S/i +V2 





Resolução: 

(a) Os valores próprios de A são dados por 

(3-A)(l-A) + l = <^= 



A 2 -4A + 4 = 
A = 2 . 



Os vectores próprios associados ao valor próprio 2 são dados por 



1 1 
-1 -1 



Como quaisquer dois vectores próprios são linearmente dependentes, escolhe-se um 
vector próprio, por exemplo, 



v = 



e procura-se um vector próprio generalizado, w 

(A - XI)w = v <==> 



1 1 
-1 -1 



a 

b 



a + b= 1 . 



Por exemplo, 



w = 



é vector próprio generalizado. Portanto, 

„At _ 



1 
-1 1 



e 2í te 2t 
e 2í 



1 
1 1 



s- 1 



(l + í)e 



2í 



te 



-te 



2t 



(l-í)e 
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(b) A solução do problema de valor inicial é dada por: 

y(t) = e At y(0) 



(l + í)e 



■2t 



te 



2í 



2í 



-te 

e 2í 
„2í 



(1 - í)e 2í 



(c) £sta equação diferencial pode ser escrita matricialmente na seguinte forma: 



Vi 
ih 



= A 



Vi 

1)2 



+ b(t) 



onde 



b(t) = 



J2t 



O 



A sua solução geral é dada pela fórmula de variação das constantes: 
y(t) = e At c+ í e A{t - s h(s)ds 



(l + í)e 



2í 



te 



+ 



2í 



(l-í)e 



Cl 
C2 



(1 + í - s)e 2 ( í " s ) (í - s)e 2 (*- s ) 

-(í-s)e 2 ( í " s ) (l-í + s)e 2 ( í - s ) 



= e 



2/ 



Ci(l + í) + C2Í 

-cií + c 2 (l-í) 



+ e 



2í 



,2í 



Cl(l + í) + C 2 í + í+ Ç 

-cií + c 2 (l-í) - f 



1 + í-s 
-í + s 



Vi e 



onde ci,C2 £ 



,2.s 



c/s 



ds 



D 



(6) 



(7) 



Considere a matriz 


5 1 " 
-1 3 




(a) Calcule e At . 


(b) Determine a solução do seguinte problema de valor inicial 




m 
. 2/2 _ 


= A 


Vi 
. 2/2 . 


+ 





com 


" yi(o) " 
. y 2 (o) _ 


= 


" " 






Considere a matriz 

A = 


2 2 " 

-2 7 




(a) Calcule e At . 


(b) Determine a solução do seguinte problema de valor inicial 




m 
. 2/2 . 


= A 


. 2/2 . 


+ 


r e 3í 



com 


" yi(o) " 
. y 2 (o) _ 


= 


" 2 " 
1 





10 
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(8) 



Considere a equação diferencial 



2/1 




4 


5 






2/1 


2/2 




-2 


-2 






2/2 


2/3 


= 






-2 




2/3 


2/4 








1 1 




2/4 


2/5 








1 1 




2/5 



(*) 



onde as entradas omitidas na matriz são zeros. 

(a) Determine a solução de (■*•) com condição inicial 

2/i (0) = 2/ 2 (0) = , 2/ 3 (0) = 2/4(0) = -y 5 (0) = 1 . 

(b) Determine a solução de (*) com condição inicial 

2/i (0) = 1 , 2/ 2 (0) = ^(0) = 2/4(0) = sfe (0) = . 

(c) Determine o conjunto de todas as condições iniciais, yo £ K 5 , tais que as 
correspondentes soluções do problema de valor inicial 

(*) 

2/(0) = 2/o 

são limitadas. 



Resolução: Seja A a matriz dos coeficientes. A solução de um problema de valor inicial 

equação (*) com y(t ) = y 

é 

y{t) = e At y Vi e R . 

Para calcular a exponencial da matriz At aproveitam-se os cálculos do exercício 3 alíneas 
(b) e (d): 

Se Ai = 



4 5 
-2 -2 



e Alí = e* 



, então 

cos í + 3 sin t 5 sin t 

—2 sin í cos t — 3 sin í 



Se A 2 = 



1 1 
1 1 



então 



,A 2 í 



e 2f + l e 2f -l 
2 2 

e 2t -l e 2f + l 
2 2 



Como 



tem-se que 





'Ai 






A = 




-2 










^2 J 



„Aí 



r e A lt 








e-' M 








e A 2 t J 



(a) A solução de (•) com condição inicial 

l/i(0) = 2/ 2 (0) = , j/ 3 (0) = 2/4(0) = -2/ 5 (0) = 1 
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11 



At 



" 









1 







e" 2í 


1 




1 


-1 




-1 



y(t) = e 



(b) A solução de (*) com condição inicial 

yi (0) = 1 , ifc(0) = 2/ 3 (0) = 2/ 4 (0) = 2/ 5 (0) = O 



Aí 



1 




cos t + 3 sin í 







—2 sin í 





= e l 






















y(t) = e 



(c) A solução de um problema de valor inicial 



equação (*) com y(0) = 



(12 

04 
«5 



!/(*) = 



e* cos í + 3e* sin í 5e* sin í 

— 2e*siní e* cos í — 3e* sin í 



-2í 



e 2t + l 
e 2t -l 



2 
e 2Í + l 



«1 
«2 
«3 
04 
a 5 



/A exponencial e envolve as seguintes funções elementares 



e cosi 



e siní 



-2í 



1 



„2í 



y4 única função limitada desta lista é a constante 1. Para que uma solução seja 
limitada, as condições iniciais, 2/(0), devem ser tais que todas as outras funções não 
apareçam. Logo, terá que ser 



0,1 = 0-2 = 0,3 = 



a 4 



-ã5 ■ 



O conjunto de todas as condições iniciais tais que as correspondentes soluções do 
problema de valor inicial são limitadas é 



' 


" 


> 









< 2/(0) = 



a 


: aeR > 


>. 


—a 


V 



D 



12 
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(9) 



Considere o sistema de equações diferenciais 








f m = - 


(siní)j/i + e cosí 








(**) 










[ m = e 


5 iní-cosí yi + ( C0S í)y 2 + 3íe siní . 


(a) 


Determine 


a solução geral do siste 


ma homogéneo 

-(siní)yi 


associado: 






[ V2 = 


3 siní-cosí yi + ( COS í)y 2 . 


(b) 


Determine 


a solução geral de (**) 






(c) 


Determine 


a solução de (**) com 


condição inicial: 








2/1(0) = 


1 e 3, 2 (0) = 


= . 



Resolução: 

(a) O sistema homogéneo pode ser resolvido em duas etapas. 
Primeiro resolve-se a primeira equação 

j/l = -(siní)yi . 



Para Vl + 0, 



— = — sin t 

Vi 



/ f- dt = - J sin t dt + 



^=> j ^ dyi = cos t + c 

<í=> In 1 2/1 1 = cosi + c 

^^ |j/i(f)| = A;ie cosí onde ki > 

^^ yi(t) = /eie cosí onde jfci ^ . 

Quando y\ se anula, tem-se que j/i(í) = 0, Vi e R, também é solução. Logo, a 
solução geral da primeira equação é 

yi(t) = kie cost , VíeR, onde fa G R . 

De seguida, substitui-se a expressão geral para y\ na segunda equação e resolve-se 
para obter y 2 : 



y-2 



s'mt—cost 
Q siní 



fcie cosí ) + (cost)y 2 



= k\e a + (cosi) y 2 . 

-a(t) 

Esta equação linear admite o factor de integração 

e I a(t)dt _ g-siní 



Mutipl içada por e smí , a equação fica 

e~ smt y 2 — (cosí)e 

, ; d 



e~ sint y 2 = ht+k 2 



Vi = ki 



V2(t) = (kit + k 2 )e 



sin t 
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A solução geral do sistema homogéneo é: 

' yi (í) = fcie cosí 

Vi eK 
y 2 (*) = (M+fc 2 )e siní , 

onde ki, k 2 £ R. 

(b) A solução geral de (**) pode ser obtida a partir da solução geral do sistema homogéneo 
pela fórmula de variação das constantes: 



yi(*) 

2/2 (í) 



= y(t) 



Cl 

C2 



+ / YWis)- 1 



3se sins 



ds 



onde Y(t) é uma solução matricial fundamental do sistema homogéneo. Obtém-se 
uma tal matriz Y(t) tomando para colunas soluções linearmente independentes do 
sistema homogéneo, como por exemplo 



Y(t) = 



e cosí 

sint j.„sint 



onde se fixou ki = 0, k 2 = 1 para a primeira coluna e ki = 1, k 2 = para a segunda. 
Com esta escolha, a inversa de Y(s) é 



Y(s)- 1 = 



-1 



3 coss+sm s 



Assim, a solução geral de (**) é 



2/1 (*) 
2/2 (*) 







O cosí 



pSÍní j-psin í 



Cl 
C 2 



+ 



e cosí 



-1 



c 2 e cosí 

Cl gSiní _|_ C2 ^ e ; 



sin t 



o e 



+ 



coss+sin s 



—e"' 

e cosí 

3 siní fpSin t 



2s 
1 



3se sin S 



(Is 



(Is 



c 2 e 



cosi 



Ci gSini _|_ c ^ e 



siní 



+ 



e cosí 

-jsiní fpSin t 



(c 2 + t)e cost 
(ci + c 2 t + 2í 2 )e siní 



Vt£ 



onde ci , C2 G 



(c) /A solução deste problema de valor inicial pode ser obtida a partir da solução geral do 
sistema homogéneo pela fórmula de variação das constantes: 



2/1 (*) 
2/2 (*) 



= r^mo)- 1 



2/1 (0) 
2/2(0) 



F / Y(t)Y(s) 

'o 



3se sins 



ds , 



14 
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onde Y(t) é uma solução matricial fundamental do sistema homogéneo, como por 
exemplo a utilizada na alínea anterior. Fazendo os cálculos, 



vi(t) 

2/2 (*) 



O e cosí 

D sin t y-psin í 



p — 1 „cosí 

e -l íe siní 



+ 



e 

1 

íe cosí 
2í 2 e siní 



té 



cosi 



2fe 



2^sin í 



(e _1 + í)e COBt 

(e- 1 í + 2í 2 )e siní 



Víe 



D 



Resolva o seguinte problema de valor inicial: 



(10) 



dy 
dt 



2 1 
2 
1 3 



e 2í 


„2í 



com 



2/(0) 



Resolução: Começa-se por resolver a segunda equação escalar que dá y 2 : 



dm 

dt 



= 2y 2 



com 



2/2(0) = . 



A solução é 

y 2 (t) = , Vi e R . 

Substituindo y 2 , resolve-se agora a terceira equação escalar que dá y%: 



dm 

dt 



3í/3 + e 



2í 



com 



2/3(0) = . 



A solução é 



VÁt)= / e^-^e 2 'ds = e 3i -e 



It 



Vte 



Finalmente, substitui-se y% e resolve-se a primeira equação escalar que dá y\\ 

C ^- = 2y 1 + e 3t com yi(0) = . 

A solução é 

Ví el . 



2/1 (*) = / e^-^e Ss ds = e st -e 2t , 

Jo 

O problema de valor inicial dado tem a seguinte solução: 



Ví e 



2/1 




" e 3í _ e 2í " 


2/2 
2/3 _ 






_ e 3t - e 2t _ 



D 

Comentário: Em alternativa, pode-se calcular a exponencial da matriz dos coeficientes 
e aplicar a fórmula de variação das constantes para sistemas de equações lineares. <) 
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(11) 



Suponha que as funções 






é + e 2í 




e í _|_ g3í 




gí _ g3í 






e 2t 


; 


e 3í 


, e 


-e 3t 











e 3í 




_ e 3t 




são três soluções y(t) da equação y 
Justifique. 


= Ay. Determine os valores próprios de A. 



Resolução: Pela fórmula de variação das constantes, as soluções de sistemas lineares de 
equações de primeira ordem homogéneas, y = Ay, são combinações lineares de exponen- 
ciais multiplicadas por potências de t da forma t k e xt , onde A é um valor próprio complexo 
da matriz dos coeficientes, A. 

As três funções dadas envolvem as exponenciais é ' , e 2t e e 3í e o sistema de equações é 
3 por 3. 

Logo, os valores próprios de A são 1, 2 e 3. D 

Equações de Ordem Superior à Primeira 



(12) 



Considere a eqi 


ação 


diferencial escalar 














y {2) + y 


= cosi . 


(*) 




(a) 


Determine 


a sol 


ução geral da equação homogénea 


associada 


a(*). 


(b) 


Determine 


uma 


solução particular de 


(*)■ 






(c) 


Determine 


a sol 


ução geral de (*). 









Resolução: 

(a) A equação homogénea associada a (*) é 

2/2) +í/ = ^ (D 2 + l)y=0 . 

O seu polinómio característico, 

p{\) = A 2 + 1 , 

tem as seguintes raízes: 

p(\) = ^^ A = ±i . 

A solução geral complexa da equação homogénea é pois 

y(t) = kie u + k 2 e~ H Vi e R 

onde ki,k2 6 C. 

A solução geral (real) da equação homogénea é 

y(t) = ci cos t + C2 sin t Vi £ R 

onde ci,C2 £ R. 

Comentário: A solução geral real obtém-se extraindo as partes real e imaginária da 

solução geral complexa, usando a fórmula de Euler: 

e lt = cos t + i sin í . 

Reorganiza-se e rebaptiza-se as constantes para simplificar a expressão final. Convém 
fazer o exercício de passagem da solução geral complexa para a solução geral real. 
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(b) Adopta-se o método dos coeficientes indeterminados: 

A função cosi é aniquilada pelo operador diferencial D 2 + 1. Se y(t) for solução 
particular de (*), i.e, 

(D 2 + l)y = cost , 
então, aplicando D 2 + 1 a ambos os membros, fica 

(D 2 + 1)(D 2 + l)i/ = (D 2 + 1) cosi = . 

Vai-se procurar uma solução particular de (*) entre a solução geral da equação ho- 
mogénea 

(D 2 + l) 2 y=0 

a qual é dada por 

y(t) = c\ cos t + ci sin t +C3Í cos t + c^t sin í . 

Os do/s primeiros termos, sobre a chaveta, são solução da equação homogénea asso- 
ciada a (*), /og"o não adiantam na pesquisa de uma solução particular da equação não- 
homogénea (*). Para encontrar uma solução particular, substitui-se y(t) = c%t cos t + 
C4ÍSÍ11Í em (*) e determina-se os coeficientes C3 e C4: 

(D 2 + 1) (cst cos í + C4Í sin í) = cos t 
<í=> — 2c3 sin í + 2c4 cos í = cos t 
<í=> C3 = e C4 = | . 

Conclui-se que uma solução particular de (•) é, por exemplo, y(t) = |ísiní para 
qualquer íeR. 

(c) Como se irara de uma equação linear, a solução geral de (*) é da forma 

solução particular \ J solução geral da equação 
de (*) J \ homogénea associada a (•) 

,4ss;m, a solução geral de (•) é 



2/(í) = ií sin í + ci cos í + C2 sin í Vi £ 



onde ci,C2 £ 



D 



(13) 



Determine a soli 


ção gera 


Id 


e cada 


uma das seguintes equações diferenciais escalares. 


(a) 








y (3) + y = , 


(b) 








y (3) + y = e* , 


(c) 








y( 3 ) + y = te* , 


(d) 








y (3) + y = 1 , 


(e) 








y (3) + y = 1 + cosi , 


(f) 








y (3) + y = e 2t cost . 
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Resolução: 

(a) Esta equação pode ser escrita 

(D 3 + D)y = 0. 

seu polinómio característico, 

p(X) = A 3 + A , 

tem as raízes 

p(\) = o ^^ A = ou A = ±i . 

A solução geral complexa da equação homogénea é pois 

y{t) = k + kie a + k 2 e~ it , Vi G R , 

onde ko,ki,k.2 G C. 

/4 solução geral (real) da equação homogénea é 

y(t) = co + ci cos í + C2 sin í , Vi G R , 

onde co,ci,C2 G R. 

Comentário: A solução geral real obtém-se extraindo as partes real e imaginária da 

solução geral complexa, usando a fórmula de Euler: 

e lt = cos t + i sin t . 

Reorganiza-se e rebaptiza-se as constantes para simplificar a expressão final. Recomenda- 
se o exercício de passagem da solução geral complexa para a solução geral real. <^ 



(b) Como se trata de uma equação linear, a solução geral é da forma 

solução geral da equação 



{ solução particular } + ■> , . , 

1 homogénea associada 

Para determinar uma solução particular, adopta-se o método dos coeficientes inde- 
terminados. A função e* é aniquilada pelo operador diferencial D — 1. Se y(t) for 
solução particular da equação dada, i.e, 

D(D 2 + l)y = ê , 

então, aplicando D — 1 a ambos os membros, fica 

(D - l)D{D 2 + 1)í/ ={D- iy = . 

Assim, vai-se procurar uma solução particular entre a solução geral da equação ho- 
mogénea 

(D-l)D(D 2 + l)y = 

a qual é dada por 

y(t) = co + ci cos t + C2 sin t +C3e* . 

Os três primeiros termos, sobre a chaveta, são solução da equação homogénea asso- 
ciada, logo não adiantam na pesquisa de uma solução particular da equação não- 
homogénea. Para encontrar uma solução particular, substitui-se y(t) = c^e 1 na 
equação e determina-se o coeficiente c^: 

D{D 2 + l)(c 3 e*) = é 
<í=> 2c3e* = é 
^^ c 3 = \ ■ 
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Uma solução particular é, por exemplo, y(t) = \ê para qualquer íeR. Conclui-se 
que a solução geral desta equação é 

y(t) = co + ci cos t + C2 sin t + -é , Vi £ M. . 

(c) Como se trata de uma equação linear, a solução geral de é da forma 

solução geral da equação 



{ solução particular } + , , . , 

1 homogénea associada 

Para determinar uma solução particular, adopta-se o método dos coeficientes inde- 
terminados. A função te t é aniquilada pelo operador diferencial (D — l) 2 . Se y(t) for 
solução particular da equação dada, i.e, 

D(D 2 + l)y = te 1 , 

então, aplicando (D — l) 2 a ambos os membros, fica 

(D - 1) 2 D(D 2 + 1)í/ ={D- 1)V = . 

Assim, vai-se procurar uma solução particular entre a solução geral da equação ho- 
mogénea 

(D-l) 2 D(D 2 + l)y = 

a qual é dada por 

y(t) = co + ci cos t + C2 sin t +cze l + C4Íe í . 

Os três primeiros termos, sobre a chaveta, são solução da equação homogénea as- 
sociada, logo não adiantam na pesquisa de uma solução particular da equação não- 
homogénea. Para encontrar uma solução particular, substitui-se y(t) = c^é" + c^e* 
na equação e determina-se os coeficientes C3 e C4 : 

D(D 2 + l)(c 3 e* + c 4 íe í ) = tê 
<==> 2c 3 e* + 4C4C* + 2c 4 íe í = té 
^^ 2c 3 + 4c 4 = e 2c 4 = 1 
<S=> c 3 = -1 e c 4 = \ . 

Uma solução particular é, por exemplo, y{t) = —e l + \tê para qualquer t £ R. 
Conclui-se que a solução geral desta equação é 

y{t) = co + ci cos t + c 2 sin t + -e t + -té , Vi £ R . 

(d) Como se trata de uma equação linear, a solução geral de é da forma 

solução geral da equação 



{ solução particular } + , , , . , 

1 homogénea associada 

Para determinar uma solução particular, adopta-se o método dos coeficientes inde- 
terminados. A função 1 é aniquilada pelo operador diferencial D. Se y{t) for solução 
particular da equação dada, i.e, 

D(D 2 + l)y=l, 

então, aplicando D a ambos os membros, fica 

D 2 (D 2 + l)y = D1 = . 

Assim, vai-se procurar uma solução particular entre a solução geral da equação ho- 
mogénea 

D 2 (D 2 + l)y = 
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a qual é dada por 

y(t) = co + ci cos t + ci sin t +C3Í . 

Os ires primeiros termos, sobre a chaveta, são solução da equação homogénea asso- 
ciada, logo não adiantam na pesquisa de uma solução particular da equação não- 
homogénea. Para encontrar uma solução particular, substitui-se y(t) = c%t na 
equação e determina-se o coeficiente C3: 

D(D 2 + 1)(C3Í) = 1 <«=> es = 1 ■ 

L/ma solução particular é, por exemplo, y(t) = t para qualquer íel. Conclui-se que 
a solução geral desta equação é 

y(t) = co + ci cos t + ci sin í + t , Vi e R . 

(e) Como se traía de uma equação linear, a solução geral de é da forma 

solução geral da equação 



{ solução particular } + , , , . , 

1 homogénea associada 

Para determinar uma solução particular, adopta-se o método dos coeficientes indeter- 
minados. A função 1 + cos t é aniquilada pelo operador diferencial D{D 2 + 1), porque 
D aniquila 1 e D 2 + 1 aniquila cosi. Se y(t) for solução particular da equação dada, 
i.e, 

D(D 2 + l)y= 1 + cosí , 

então, aplicando D(D 2 + 1) a ambos os membros, fica 

D 2 {D 2 + l) 2 y = D{D 2 + 1)(1 + cosi) = . 

Assim, vai-se procurar uma solução particular entre a solução geral da equação ho- 
mogénea 

D 2 (D 2 + l) 2 y = 
a qual é dada por 

y(t) = co + ci cos t + C2 sin t +C3Í + C4Í cos t + c$t sin t . 

Os três primeiros termos, sobre a chaveta, são solução da equação homogénea as- 
sociada, logo não adiantam na pesquisa de uma solução particular da equação não- 
homogénea. Para encontrar uma solução particular, substitui-se y(t) = C3Í+C4ÍCOSÍ+ 
csísiní na equação e determina-se os coeficientes C3, C4 e C5: 

D(D 2 + l)(c3Í + C4ÍC0SÍ + csísiní) = 1 + cosi 
<í=> C3 — 2c4 cos t — 2c5 sin í = 1 + cos t 
<í=> C3 = 1 e C4 = — j ec5 = 0. 

Uma solução particular é, por exemplo, y(t) = t — | cos í para qualquer t 6 R. 
Conclui-se que a solução geral desta equação é 

y(t) = co + ci cos í + C2 sin í + t — — cos í , Vi £ R . 

(f ) Basía encontrar uma solução particular para a equação. A função e 2t cos t é aniquilada 
pelo operador diferencial (D — 2) 2 + 1. ,4ss/m, uma solução particular da equação 
será uma solução da equação homogénea 

D(D 2 + l)((D-2) 2 + l)y = 

A solução geral desta equação é dada por 

y{t) = co + ci cos t + C2 sin t + c%e 2t cos t + c^e 21 sin t 
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Os três primeiros termos são soluções da equação homogénea inicial por isso podemos 
fazer cq = c\ = C2 = 0. Tendo em conta que 



„2í. 



(D 2 + 1) (eH cos t) = 4e 2i cos t - 4e 2í sin t 



e que 



{D 2 + 1) (e 2í sin st) = 4e 2í cos t + 4e 2í sin t , 
substituindo c 3 e 2í cosí + C4e 2í siní na equação, obtém-se 

(D 2 + l)D(c 3 e 2t cos í + c 4 e 2í sin í) = e 2 t cos í 

(D 2 + 1) ((2c 3 + c 4 )e 2í cos í + (2c 4 - c 3 )e 2í sin í) = e 2 t cos t 

(4(2c 3 + c 4 ) + 4(2c 4 - c 3 ))e 2í cos í + (4(2c 4 - c 3 ) - 4(2c 3 + c 4 ))e 2í sin í = e 2 t cos í 

4c 3 + 12c 4 = 1 
4c 4 - 12c 3 = 

° 3 = ¥ 
c 4 = íõ 



Logo, 



y(í) = -— e 2í cos t + —e 2t sin t 
yy ' 40 40 



é t7A7ia solução particular. Conclui-se que a solução geral é 
y(t) = cot + ci cos t + C2 sin í + 



•L 9+ *J 9+ 

— e cosí+ — e smí , 

40 40 



Víe 



D 



(14) 



Determine a 


solução gera 


1 d 


e cada uma das seguintes equações diferenciais. 


(a) 






y (3) + y (2) + y +y = 0; 


(b) 






y(3) _|_ y (2) _|_ y _|_ y = l _|_ gí gin ^ . 


(c) 






y (2) -y- 2 y = 2e 2t - 2 ; 


(d) 






y (3) _ 27ry ( 2 ) + 7T 2 y = 27T 3 Í ; 


(e) 






y ( 3 ) _ 2y (2) + 2y = e"*cosí 



(15) 



Determine a solução do seguinte problema de valor inicial: 



y(3) + y 



y(0) = 1 , y(0) = 1 , ^(0) = | 



Resolução: Pe/a alínea (b) do exercício 13, a solução geral da equação diferencial é 
y(t) = co + ci cos t + C2 sin í + -e* Vi 6 R . 
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As derivadas desta solução são 

y(t) = — ci siní + c 2 cosi + \ê 

y^ 2 \t) = —ci cosi — c 2 siní + \ê . 
A condição inicial impõe que 

1 = 2/(0) = co + ci + \ 

1 = 2/(0) = c 2 + \ 

\ = 2/ {2) (0) = -ci + i 

donde se conclui que 

1 1 

co = 2 > ci = e c 2 = - . 

Logo, a solução deste problema de valor inicial é 

V(t)=- + -smt+-e t , VteM. 



D 



(16) 



Determine a solução da equação linear escalar 




yW + 2 y( 2 ) + y = b(t) 




que verifica as condições iniciais y(0) = y(0) = 0, y^^O) = 


= 1, quando 


(a) b(t) = 0, Vi G R; 




(b) 6(í) = í, Vi G K; 




(c) 6(t) = e*. Vi G K. 





Resolução: 

(a) /\ equação homogénea pode ser escrita 

(D 3 + 2D 2 + £>)y = ou seja D(D + l) 2 y = 
cu/a solução geral é 

y(t) = co + cie _í + c 2 íe"* Vi e R 

onde co, ci, C2 G R. /4s derivadas desta solução são: 

y(t) = -cie' 1 + c 2 e _í - c 2 íe _í 
?/ (2) (í) = cie _í -2c 2 e _í + c 2 íe _í . 

A/o va/or inicial tem-se 

2/(0) = co + ci 

2/(0) = -ci + c 2 

2/( 2 )(0) = ci-2c 2 . 

/\ condição inicial impõe que 



CO + Cl = 


= 


[ Co = 


1 


-ci + c 2 = 


= <^= 


> Cl = 


-1 


ci - 2c 2 = 


= 1 


l -2 = 


-1 



Z-ogo, a solução do problema é 

y(t) = 1 - e"* - íe _í Vi G 



22 AMIV - FICHA SUPLEMENTAR 4 

(b) A solução geral de uma equação linear não homogénea pode ser obtida somando 
uma solução particular ã solução geral da equação homogénea associada. Para obter 
uma solução particular da equação com b(t) = t, aplica-se o método dos coeficientes 
indeterminados: 

-té solução de D 2 y = 0; 

- se y é solução de (D 3 + 2D 2 + D)y = t, então 

D 2 (D 3 + 2D 2 + D)y = D 2 t=0 ; 

- procura-se uma solução particular da equação dada entre a solução geral da 
equação homogénea, 

D 2 (D 3 + 2D 2 + D)y = ou seja D 3 (D + l) 2 y = , 

que é 

y(t) = co + cie - * + C2Íe~*+c 3 í + c 4 í 2 ; 

- os termos sobre a chaveta constituem a solução geral da equação homogénea, 
logo não adiantam na busca de uma soluçãoparticular da equação não ho- 
mogénea; 

- toma-se como candidata para solução particular uma função da forma 

y{t) = c 3 t + c 4 í 2 
a qual tem as seguintes derivadas 

y(t) = C3 + 2c4Í 

y {2 \t) = 2 C4 
y {3 \t) = 0; 

- os coeficientes cy, e c 4 determinam-se substituindo na equação: 

y (3) + 2y W + y = t 

<==> + 4c 4 + c 3 + 2c 4 í = t 
<í=> 4c 4 + c 3 = e 2c 4 = 1 
<í=> c 3 = -2 e c 4 = \ . 

Obteve-se a solução particular 

y(t) = -2t+-t 2 VíeM. 

A solução geral da equação diferencial dada é 

y(t) = -2í + \t 2 + co + cie _í + c 2 íe _í Vi e R . 

sol. particular sol. geral da eq. hom. 

onde co, ci, C2 e R. As derivadas desta solução são: 

y{t) = -2 + í-cie- í + c 2 e- í -c 2 te- í 
y {2) (t) = l + c 1 e- t -2c 2 e- t + c 2 te- t . 

No valor inicial tem-se 

2/(0) = co + ci 
2/(0) = -2-ci + c 2 
2/( 2 ) (0) = l + ci-2c 2 . 

/\ condição inicial impõe que 

co + ci = ( co = 4 

-2-ci + C2 = ^^ < ci = -4 
1 + ci - 2c2 = 1 | c 2 = -2 
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Logo, a solução do problema é 

y(t) = -2í+ -í 2 + 4-4e" í -2íe" í Vi e R . 

(c) 4 solução geral de uma equação linear não homogénea pode ser obtida somando uma 
solução particular ã solução geral da equação homogénea associada. Para obter uma 
solução particular da equação com b(t) = ê , aplica-se o método dos coeficientes 
indeterminados: 

- é é solução de (D — l)y = 0; 

- se y é solução de (D 3 + 2D 2 + D)y = é , então 

(D - 1)(D 3 + 2D 2 + D)y = (D - l)e* = ; 

- procura-se uma solução particular da equação dada entre a solução geral da 
equação homogénea, 

(D-1)(D 3 + 2D 2 + D)y = 

ou seja 

(D-l)D(D+l) 2 y = 0, 

que é 

y(t) = co + cie _í + c 2 íe _í +c 3 e t ; 

- os termos sobre a chaveta constituem a solução geral da equação homogénea, 
logo não adiantam na busca de uma soluçãoparticular da equação não ho- 
mogénea; 

- toma-se como candidata para solução particular uma função da forma 

y(t) = c 3 e* 
que tem as seguintes derivadas 

y(t) = y {2) (t) = y^(t) = c 3 e t ; 

- o coeficiente c 3 determina-se substituindo na equação: 

y (3) + 2y (2) + y=e t 



(ca + 2c 3 + c 3 )é = é 
Obteve-se a solução particular 



C3=l 



y(t) = \ê , Vi e R . 

A solução geral da equação diferencial dada é 

y(t)= \é + cp + cie~ * + c 2 te~* Vi e 

sol. part. sol. geral da eq. hom. 

onde co, ci, C2 e R. /4s derivadas desta solução são: 

y(t) = \e l - c\e~ l + c 2 e _í - c 2 íe _í 
2/ (2) (*) = je* + cie _í - 2c 2 e _í + c 2 íe _í . 

No valor inicial tem-se 

2/(0) = ^ + c + ci 

y(o) = j - ci + c 2 

2/ (2) (0) = i + ci - 2c 2 . 
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A condição inicial impõe que 



( 1 

4 
1 



+ co + ci = 
- ci + c 2 = 



v \ + ci-2c 2 = 1 
Logo, a solução do problema é 



-te" 



Cl 

Vi e 



= 



D 



(17) 



Determine a solução do seguinte problema de valor inicial 

2y(2) _|_ y = _4 _|_ 2í 



.y 



(3) 



2/(0) = 1 , 2/(0) = 0, y^(0) = 2 



(18) 



Considere 


a equação 


diferencial esca 


lar 




















y W 


+ y (3) 


+ y (2) 


+ y = 


. 








(a) 


Determine a sua 


solução geral. 
















(b) 


Determine para 


que condições 


iniciais 


em t 


= é 


que os 


problemas 


de 


valor 




inicia 


correspon 


dentes têm soli 


jção convergente quando t 


— > +00. 







Resolução: 

(a) O polinómio característico da equação é 

p(X) = A 4 + A 3 + A 2 + A 

cuja factorização em monómios é 

p(A) = A(A + l)(A-i)(A + i) . 

Comentário: Para chegar a esta factorização, repara-se na raiz e adivinha-se a 
raiz-1. 

A solução geral (real) da equação é 

y(t) = ci + C2e _í + C3 cos t + C4 sin t Vi e R 

onc/e ci, C2, C3, C4 G R. 

(b) Para que uma solução seja convergente quando t —> +00, ela não pode envolver 
as funções cosi nem siní. Procuram-se então as condições iniciais em t = que 
implicam que es e 04 sejam 0. Uma vez que 

y(t) = ci + C2e _í + C3 cosi + C4siní 

y(t) = — C2e~* — C3 sin í + C4 cosi 

V (*) = C2e _í — C3 cosi — C4 sin í 

y( 3 \t) = — C2&~ t + C3 sin t — C4 cosi 

os valores em t = são 

í/(0) = Cl + C2 + C3 

y(0) = - c 2 + c 4 

y (2) (o) = C2-C3 

y (3) (o) = -C2-C 4 
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donde sai 


que 












Cl = 
C2 = 
C3 = 

C 4 = 


?/(0) + 2/(0) -2/ 2 )(0) + y( 3 ) (0) 
-i(2/(0) + í/( 3 )(0)) 
-^(2/(0) + í/ {3) (0))-|/ (2) (0) 
i(2/(0)-2/ 3 )(0))- 


Para que 
terá que > 


a solução do problema de valor inicial seja convergente quando t 

ser 






= 
= 


-è(y(o) + y (3) (o))- 

i(2/(0)-í/( 3 )(0)) 


-2/ (2) (0) 



+oo, 



ou seja, 

2/(0) = -2/ (2) (0) = 2/ {3) (0) . 
O conjunto de todas as condições iniciais tais que as correspondentes soluções do 
problema de valor inicial são convergentes quando t —> +oo é 

{(2/(0), 2/(0), 2/ (2) (0),2/ (3) (0)) = (a, 6, -6, 6) : a, b e R} . 

D 



